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Abstract 

We investigate how entanglement spreads along small Bose-Hubbard chains, with only the first 
well initially occupied by a mesoscopic number of atoms, as the number of sites increases. For two- 
and three-well chains in the non-interacting case, we are able to obtain analytical solutions and 
show that the presence of entanglement depends on having a sub-Poissonian state of the atoms 
in the first well. In these cases, the correlations we calculate are completely periodic. Restoring 
the collisional interactions or moving to a four-well chain necessitates a numerical treatment, for 
which we use the fully quantum positive-P representation. We examine two different correlations 
and find that adding collisional interactions destroys the periodicity of the correlations and causes 
them to degrade with time. This happens well before there is a noticeable effect on the periodicity 
of the solutions for the number of atoms in each well. 

PACS numbers: 03.67.Mn,03.65.Ud,03.75.Gg,67.85.Hj 


1 


I. INTRODUCTION 


The investigation of non-local correlations such as entanglement and Einstein-Podolsky- 
Rosen (EPR) steering [l-3| necessarily belongs to the held of continuous-variable entangle¬ 
ment, which is a very active area of research 5|. Many criteria having been developed 
to signify the presence of inseparability, entanglement and the potential for EPR steering, 
especially in bipartite systems. The measures that we use in this work were developed from 
work by Hillery and Zubairy jo] which was based on the Cauchy-Schwarz inequality. The 
original inequality covered bipartite inseparability and was expanded on by Cavalcanti et 


. Q 


al. (tI to cover multoartite entanglement, steering, and violations of Bell inequalities. As 
shown by He et al. |8|, the Hillery and Zubairy criteria are well suited to number conserving 
processes such as those of interest here, but as shown by Olsen jsl, they will sometimes miss 
actually existing entanglement in quantum optical systems. In previous work on a different 
three-well system, we have found these to be the most useful criteria for the detection of 
entanglement in Bose-Hubbard systems fli Q 

Multi-mode entanglement in Bose-Einstein condensates (BEC) has been predicted and 


examined in the processes of molecular dissociation 
tice 


13 


15| . and in a two-well Bose-Hubbard model 


A 


a 


our-wave mixing in an optical lat- 


of the modes is produced by the tunne 
and pulsed tunneling configurations 


16|. In the latter case the separation 

r nr 

ing between wells, in both the continuous [8|, Il7|, [18 1 
20j . The spreading of correlations and entanglement 


after a sudden parameter change, investigated by Lauchli and Kollath 2l|], and the deve 


opment of entanglement between the two ends of a chain, analysed by Reslen and Bose 22| . 
are particularly relevant to the present work. 

The quantum correlations necessary to detect inseparability, entanglement and EPR 
steering can in principle be measured using the interaction with light 23|], or by homodyning 


with other atomic modes 


24| . We note here that the entanglement we are examining is a 


collective property between atomic modes which are spatially separated, and is not between 
individual atoms [l^ 16 1. We also note that tunneling of atoms between wells will not, by 
itself, lead to ent ang lement, with either a nonlinear interaction or initial non-classical states 
being necessary [25 1. 

What we investigate in this article is the spreading of entanglement and the potential 
for EPR steering along small Bose-Hubbard chains, beginning with a mesoscopic population 
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of atoms in one end well. In the case of non-interacting atoms and two- or three-well 
chains, we will obtain analytical solutions of the Heisenberg equations of motion for the 
atomic annihilation and creation operators. In these cases, we hnd that the presence of 
entanglement or EPR steering will depend on the presence of a non-classical quantum state 
of the initial atoms. For four wells, and in all cases with atomic interactions, we will use 
numerical stochastic integration in the positive-P representation 2^, which proved to be 
stable over the evolution times we examine here. 


II. PHYSICAL MODEL AND HAMILTONIANS 


We follow the approach to the Bose-Hubbard model I2h. l28l| taken by Milburn et al. 29l |. 


also generalising this to three [30|, 


3l| and four wells [32|. All our systems are in a linear 


conhguration with only the hrst well initially occupied. Where possible, we analytically solve 
the Heisenberg equations of motion and in other cases we use the fully quantum positive-P 


representation 


2q . which contains the full dynamics due to the Hamiltonians. We consider 


these to be the most suitable approaches here because they are both exact, allow for an 
easy representation of mesoscopic numbers of atoms, can be used to calculate quantum 
correlations, and can simulate different quantum initial states 33|]. Just as importantly, 
both calculations scale linearly with the number of sites and can in principle deal with any 
number of atoms. One disadvantage of the positive-P representation is that the integration 
can show a tendency to diverge at short times for high collisional nonlinearities 3^ . As long 
as the procedures followed to derive the Fokker-Planck equation for the positive-P function 
are valid [3^, the stochastic solutions are guaranteed to be accurate wherever the integration 


converges. With all the results shown here, we found no sign of numerical divergences. 

The systems we investigate are very simple, with the potential wells in a linear conhg¬ 
uration. Fach of these can contain a single atomic mode, which we will treat as being in 
the lowest energy level. Atoms in each of the wells can tunnel into the nearest neighbour 
potential. With the Oj as bosonic annihilation operators for atoms in mode j, J representing 
the coupling between the wells, and y as the collisional nonlinearity, we may now write our 
Hamiltonians. Following the usual procedures 29|, we hnd the two-well Hamiltonian as 

2 

7^2 = , ( 1 ) 

J = 1 
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the three-well Hamiltonian as 


■Hs = + 02^1 + «3®2 + 02^3 j , 

i=i 


( 2 ) 


and that for the four well system as 


4 

n 4 = h E Xa}j -|- hJ (^d\a 2 + a^oi -|- 03^2 -|- a^cis -|- a^os -|- dld^ . 


(3) 


i=i 


From these three Hamiltonians we can derive the equations of motion used to investigate 
the dynamics of our systems. 


III. NON-INTERACTING CASE 

For the case where the collisional interaction between the atoms is set to zero, we hnd 
that an analytical solution of the Heisenberg equations of motion of the system operators 
for the two- and three-well chains is possible. For the four-well chain, we were unsuccessful. 
The two-well system is rather simple and has been extensively analysed. We will present 
the results here for completeness and for the purposes of easy comparison with the longer 
chains. 


A. Two wells 

For the two-well system, the Heisenberg equations of motion are found as 
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This set of linear operator equations is readily solved, having the solutions 

di{t) = ai( 0 ) cos — ^2(0) sin Jt, 
a|(t) = a|( 0 ) cos -|-*02(0) sin Jt, 

02(0 = —fai(O) sin -|- 02 ( 0 ) cos Jt, 
dl{t) = ia|(0) sin Jt-|- 5 , 2 ( 0 ) cos 


(4) 


( 5 ) 
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For the populations, we find the time dependent solutions, 


(iVi(t)) = (aJ(O)ai(O)) cos^ Jt, 
{N2{t)) = (a|(0)ai(0)) sin^ Jt, 


( 6 ) 


showing that the atoms will cycle forever between the two wells, as long as any other consid¬ 


erations are absent. The inclusion of t 
eventual revival of the oscillations 
here. 


le atomic interactions is known to cause a collapse and 


3l| . but over longer timescales than we are investigating 



5i 
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B. Triple wells 

For the non-interacting three-well system, the Heisenberg equations of motion are 

ai(0) 

aUo) 

02(0) 
o5(o) 

03(0) 

4 ( 0 ) 

This set of linear operator equations is also readily solved, having the solutions 

1 i 1 

ai{t) = - (cosHt1) ai(0)-^ sin02(0)- (cos — 1) 03(0), 

2 ^/2 2 

1 z 1 

aj(f) = - (cosHt -|- 1 ) a|( 0 ) -I—^ sin 02(0) H— (cos — 1 ) 03(0), 

2 \/2 2 

—z z 

02(t) = —;= sin ai( 0 ) -|-COS 02(0) -^ sin 03(0), 

\/2 V 2 

al{t) = —;= sin fit a|^( 0 ) + cos fit 5.2(0) H—^ sin fit 03(0), 
v 2 V 2 

1 i 1 

asi^t) = - (cosHf — 1 ) 5 i( 0 )-^ sin fit 02(0)- (cos fit-|- 1 ) 03(0), 

2 1/2 2 

1 z 1 

al{t) = - (cosHf — 1) aj(0) H—^ sin 02(0)- (cos -|-1) 03(0), 

2 1/2 2 


( 7 ) 


( 8 ) 


where we have made the substitution H = \/2 J for reasons of notational elegance. In this 
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FIG. 1: (Color online) The analytical non-interacting solutions for the populations in each well as 
a function of time, for J = 1 and A^i(O) = 200, with A^2(0) = iV3(0) = 0. The interacting stochastic 
results for x = 10“^ and both coherent and Fock initial states are indistinguishable on this scale. 
The quantities plotted in this and subsequent plots are dimensionless. 

case we find the populations as 

= ^(cosQt + if (al(O)ai(O)}, 

(N 2 (t)} = ^sin^Gt(a|(0)ai(0)), 

(Nsit)) = ^ (cosGt - 1)^ (a|(0)ai(0)). (9) 

These solutions again describe fully periodic oscillations and are shown in Fig. [1] 


Four wells 


In the four well case, we were not able to obtain analytical solutions of the Heisenberg 
equations of motion. However, we can integrate the positive-P representation 2^ equations 
numerically with initial states taken from distributions using the methods developed by 
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FIG. 2: (Color online) The numerical non-interacting solutions for the populations of the four-well 
model as a function of time, for J = 1, and iVi(O) = 200 in a Fock state, with A^2(0) = N^^O) = 
Ni(0) = 0. These solutions were averaged over 2.88 x 10^ stochastic realisations of the positive-P 
equations. We see that the populations now exhibit much more complicated dynamics than for the 
three-well case shown in Fig. [TJ 


Olsen and Bradley 33j. For an initial coherent state, this is the same as integrating the 
classical, Gross-Pitaevskii equations. For an initial Fock state, while each trajectory is 
deterministic, the initial state is chosen from a distribution in the doubled phase space of 
the representation. The equations are as given in section |V] below, but without the noise 
terms. We present the solutions for the populations here, beginning again with an initial 
Fock state of 200 atoms in the hrst well. 


IV. ANALYTIC QUANTUM CORRELATIONS 

As well as the populations in each well, we also calculate any type of operator products 
that we desire, analytically in the case without interactions. We will calculate the number 
variances in each well, and entanglement and EPR steering measures adapted from an 
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inequality developed by Hillery and Zubairy. They showed that, considering two separable 
modes denoted by i and j j6|, 

\{a\aj)\^ < {a\aia]aj), ( 10 ) 

with the equality holding for coherent states. The violation of this inequality is thus an 
indication of the inseparability of, and entanglement between, the two modes. Cavalcanti 
et al. Q have extended this inequality to provide indicators of EPR steering 1-3| and Bell 
violations j^|. We now dehne the correlation function 




t\ 


at; 


( 11 ) 


for which a positive value reveals entanglement between modes i and j. We easily see that 
^ij gives a value of zero for two independent coherent states and a negative result for two 
independent Fock states. 


L Q 


Cavalcanti et al. [7i further developed the work of Hillery and Zubairy to hnd inequalities 
for which the violation denotes the possibility of EPR steering and Bell states. The EPR 
steering inequality for two modes is written as 


|(oia])|2 < {alai{a]aj +-)), 


while the Bell state inequality is written as 


Ka.a;)p < ((o?a, + + 1)), 


( 12 ) 


(13) 


J' I — ' 2 ^ ^ 

Calling on the overworked Alice and Bob, if Alice measures mode i and Bob measures mode 
i a violation of the inequality (IT^ signihes that Bob would be able to steer Alice, and 
vice versa for a swapping of the modes. These inequalities allow us to dehne a correlation 
function which signihes the presence of EPR steering when it has a value of greater than 
zero, 


Sij = {aia]){a\aj) - {a\ai{a]aj + -)), 


(14) 


“J/ ■ 2' 

and another for which a positive value signihes the presence of Bell correlations, 

1 1 

Cij = {aia]){a\aj) - {{0.% + ^){a]aj + -)). (15) 

The criteria ^ij and Sjj have been shown to detect both inseparability and asymmetric EPR 
steering in a three-well Bose-Hubbard model under the process of coherent transfer of atomic 
population (CTAP) 2, 20|, as well as bipartite entanglement in a three-mode model with 
all population initially in the central well QQ- 
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A. Twin well model 


The first class of correlations we calculate are the number variances in each well. In terms 
of the operators, these are 

V{Nj) = (OjajOjaj) — (16) 


We hnd, dropping the time argument for simplicity, 

V{N^) = V{Ni{0))cos\Jt+^{Ni{0))sm‘^2Jt, 
V{N 2 ) = V{Ni{0))sm^ Jt+^{Ni{0))sm^2Jt. 


The Hillery-Zubairy correlation is found as 


62 


sin^ 2Jt 


mo)) - v"(^i(o)) 


(17) 


(18) 


showing that entanglement will be detected by this measure whenever the initial state of 
the atoms in the hrst well is sub-Poissonian. The correlation will be maximised for an intial 
Fock state, as shown in Fig. ID 

The EPR steering correlations are solved as 


Sio = -sim 2 Jf 


S 21 = - sim 2Jt 


(iVi(O)) - R(W(0))J - - sin^ Jt(iVi(0)), 
(iVi(O)) - R(iVi(0))] - ^ cos^ Jt(iVi(0)), 


(19) 


both of which can be positive. We note that S 12 7 ^ S 21 , which leaves open the possibility of 


asymmetric steering 


37|, which we see in Fig. |3l where S 12 and S 21 are positive at different 


times for an initial Fock state. We note here that, although the steering is asymmetric 
by this measure, this does not mean this is the case for all possible measures. The Bell 
correlation is found as 


C12 


sin^ 2Jt 


(iVi(o)) - r(JVi( 0 )) 



which obviously can never be positive. 


( 20 ) 
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FIG. 3: (Color online) The analytical non-interacting solutions for ^12 and Sjj in the twin-well 
model as a function of time, for J = 1 and A^i(O) = 200 in a Fock state, with A^2(0) = = 0. 

We see that the exhibit positivity at different times, showing that asymmetric steering is 
possible by this measure. 

B. Triple well model 


In the non-interacting case and again with only the hrst well initially occupied, we hnd 

'1 


V{Ni) = —(1-f cosGt)'^ y(iVi(0))-I—sin^ Gt 
16 8 


■ sim Gt -|- (1 + cos ^Ity 


(A^i(O)), 


V{N2) = - sin^ Gt I/(iVi(0))-I-- sin^ Gt — - sin^ GtJ (iVi(O)) 

V{Ny = 4(cosGt-l)V(iVi(0)) + ^(cosGt-l)2(iVi(0)). 
ib 4 


( 21 ) 


The Hillery-Zubairy inseparability correlations ^ij are found as 


■^12 


■sin2Gt(l + cosGt)^ (^i(O)) - V{Ni{0)) 


1 


16 


sin^ fit 


(iVi(o)-y(iVi(o)) 


( 22 ) 
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FIG. 4: (Color online) The analytical non-interacting solutions for ^12 and ^13 in the triple-well 
model as a function of time, for J = 1 and A^i(O) = 200 in a Fock state, with iV 2 ( 0 ) = A^ 3 ( 0 ) = 0. 
We see that the maximum of ^12 is greater than the maximum of ^ 13 , suggesting that entanglement 
decreases with distance along the chain. 


from which we can see that both modes 1 and 2, and 1 and 3 can exhibit a time dependent 
inseparability for an initial snb-Poissonian state in the hrst well. We also see from both the 
eqnations and Fig. |4] that the maximnm valne of the fnnction is less for ^13 than for ^ 12 , 
indicating that the entanglement with the hrst well decreases with distance along the chain. 

For the EPR steering criteria we hnd 


Fi2 

^21 

Fi3 

S31 


-- (1 + cosQtf (iVi(O)) cos^ nt + ^(^(0)) 

8 L 

~ sin^ fit \{1 + cos QtY ^(^(O))-I-(iVi(O)) + 2 cos^ fit{Ni{0)) \ , 


16 

1 


■ sin^ fit 


(iVi(o)) - ^(^1(0)) - - (1 + cosntf (iVi(o)) 


— sin^nt (W(0)) - ^(^ 1 ( 0 )) - - (cosflf-l)^(Ali(0)), 
16 ° 


(23) 


none of which can be greater than zero for initial coherent states. For initial Fock states, 
we hnd that S 13 and S 31 can take on positive valnes, and we again see that the correlations 
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FIG. 5: (Color online) The analytical non-interacting solutions for the Sjj in the triple-well model 
as a function of time, for J = 1 and A^i(O) = 200 in a Fock state, with N2{0) = N^lO) = 0. We see 
that only S 13 and S 31 take on positive values. 

measured are not equivalent under an exchange of indices. These are shown in Fig. [5l 


V. STOCHASTIC METHODS 


In this case (x 7^ 0); if is not obvious how to solve the full quantum equations of motion 
analytically. We will therefore use the positive-P representation 2^, which allows for exact 
solutions of the dynamics arising from the Hamiltonians, in the limit of the average of an 
inhnite number of trajectories of the stochastic differential equations in a doubled phase- 
space. In practice we obviously cannot integrate an inhnite number of trajectories, but have 
used numbers large enough that the sampling error is within the line thicknesses of our 
plotted results. Following the standard methods [3^, the set of ltd stochastic differential 
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equations [3^ for the two-well system are found as 

dai 


= -2ixafal - iJa2 + \J-2ixal rji, 
= 2ixat‘^ai + iJa^ + \J2ixaf ^ r]2, 
^ = -2ixa^a2 - + y^-2f^r73, 

= 2ixat^a2 + iJat + ^J2ix^ Vi, 


(24) 


where the rjj are standard Gaussian noises with = 0 and Vjif)Vkit') = SjkS{t — t'). As 
always, averages of the positive-P variables represent normally ordered operator moments, 
such that, for example, —)■ (af^a™). We also note that aj = (a^)* only after 

taking averages, and it is this freedom that allows classical variables to represent quantum 
operators. 

The equations for the triple and quadruple well systems are obvious extensions of Eq. [2ll 
but we will present them here for the sake of completeness. For the three-well system we 
hnd 


dai 

dt 

daf 

dt 

doi2 

dt 

da2 

dt 

das 

dt 

da^ 

dt 


—2ixafal — iJa2 -l- \J—2ixa\ t]i, 
2ixaf‘^ai iJa^ 2ixa'l ^ 772 , 
- 2 ixa'^al - iJ (ai -F a^) + - 2 ixal 773, 
2ixa^^Q;2 + {ctt + + \j2ixd^ Vi, 

-2ixatal - iJa2 sj-2ixal 

2ixat'^<^3 + iJat + \J2ixat ^ Ve, 


(25) 
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while the equations for the four-well system are 
dai 


—2ixciia\ — iJa2 + \J —2ixGti ?7i, 
2ixot‘^ai iJa.2 + \J2ix(yt ^ V 2 , 
-2ixatal - iJ (ai 03 ) \J-2ixal r]3, 
2ixa^^Q;2 + iJ («i^ + «^) + ^4, 

-2ixoL^o?3 - iJ (02 + 0 ( 4 ) + \J-2ixoil r]^, 

dcx~^ j 

^ - 2ixa^‘^OL3 iJ {a^ + at) + \j2ixat‘^ rje, 

-2ixatal - iJa^ \J -2ixal rjj, 

2ixat^a4 + iJat + 2ixat ^ 


dt 

daf 

dt 

da2 

dt 

dat 

dt 

das 

dt 


dt 

da4 

dt 

dat 

dt 


Vs- 


(26) 


These systems of equations must be solved numerically, for which we use Matlab. This allows 
us to average over a sufficient number of trajectories in a reasonable time, in most cases less 
than 2 hours. If we wished to add wells, the growth in computational time would be linear 
in the number, and we found that 4 wells did not extend our computational resources to 
any significant extent. 


VI. NUMERICAL SOLUTIONS 

For our results in the interacting case, we have chosen a nonlinearity of x = 10“^, again 
with either a Fock or coherent state with an average of 200 atoms in the first well and all 
the others initially empty. These different quantum states are simulated using the methods 
found in Olsen and Bradley 33| . We have chosen coherent states because they are the initial 
state most often used in numerical simulations, and Fock states because we consider this 
the most natural state for atoms in an isolated well. 


A. Two wells 

For the two-well system we found evidence of both entanglement and EPR steering in the 
non-interacting case. We find that the results with a collisional nonlinearity and an initial 
Fock state are initially similar, as shown in Fig. |6l but that the correlations degrade with 
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FIG. 6: (Color online) The numerical interacting solutions for ^ 12 , S 12 and S 21 in the twin-well 
model as a function of time, for J = 1, x = 10“^ and A^i(O) = 200 in a Fock state, with A^2(0) = 0. 
We see that the correlations are no longer periodic as in the non-interacting case of Fig. [3l and 
that the two EPR steering correlations give positive signals at different times. These solutions are 
averaged over 1.9 x 10® stochastic realisations of the positive-P equations. 


time. For an initial coherent state in the hrst well, we fonnd no evidence of EPR steering, 
and only a very small positive valne of ^12 ~ 0.2 at Jt ^ 2.2. We note that onr resnlts are 
not inconsistent with those of He et al. Q|, who fonnd entanglement in the gronnd state 
of a two-well system. Onr system, beginning with only one of the wells initially occnpied, 
effectively nndergoes a sndden parameter change from J = 0 to J = 1 at f = 0, and in 
that respect has more in common with the work of Lanchli and Kollath 2l|, who stndied 
the time evolntion of entanglement following a qnench. In onr case, the sndden change in 
parameters is the tnrning on of the tnnnelling at t = 0, so that onr system is then far from 
the eqnilibrinm state of the mnlti-well systems. 
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FIG. 7: (Color online) The numerical interacting solutions for ^ 12 , S 12 and S 21 in the three- 
well model as a function of time, for J = 1, x = 10“^, and A^i(O) = 200 in a Fock state, with 
-^ 2 ( 0 ) = ^'" 3 ( 0 ) = 0. We see that the correlations degrade with time when compared to the non¬ 
interacting case of Fig. U and that only one of the EPR steering correlations gives a positive 
signal for any of the times considered. These solutions were averaged over 6.74 x 10^ stochastic 
realisations of the positive-P equations. 

B. Three wells 

With an initial Fock state in the hrst well, we hnd signals of entanglement with both ^12 
and ^ 13 , as shown in Fig. [7] and Fig. [HI These initially follow the analytical non-interacting 
solntions, bnt again become different with time. We find asymmetric EPR steering in both 
cases, with only S 21 and S 31 attaining positive valnes, while S 12 and S 13 remain negative. 

For an initial coherent state, we see no evidence of EPR steering in the three-well model, 
and only a very small signal of entanglement between wells 1 and 2, with (^12 < 2.5, and 
even that only for very short intervals of time. (^13 had a maximnm valne of zero over the 
length of time of onr investigations. 
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FIG. 8 : (Color online) The numerical interacting solutions for ^ 13 , S 13 and S 31 in the three- 
well model as a function of time, for J = 1, x = 10“^ and iVi(O) = 200 in a Fock state, with 
-^ 2 ( 0 ) = A^ 3 ( 0 ) = 0. We see that the correlations degrade with time when compared to the non¬ 
interacting case of Fig. lU and that only S 31 gives a positive EPR steering signal, and this only for 
a short time. These solutions were averaged over 6.74 x 10^ stochastic realisations of the positive-P 
equations. 

C. Four wells 


The solutions for the populations of the four-well system become more irregular, without 
the clear periodicity of the smaller chains, as was shown above in Fig. |2]for an initial Fock 
state in the non-interacting case. The solutions when we begin with either Fock or coherent 
states are indistinguishable over the time interval investigated. When we add interactions, 
the solutions begin to deviate from those without interactions around Jt = 4, with the devi¬ 
ations increasing with time. This is common with Bose-Hubbard models and has been seen 
previously in two-, three-, and four well systems 29|, 1 m|- The irregularity and lack of clear 


periodicity may be a sign of the onset of chaos, which has been predicted in longer chains, 
but with much smaller well occupation numbers j^ . 40|. Irregular behaviour has previ- 
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ously 

rates 


Deen seen in a fonr-well system in a square configuration with two different tunnelling 


31 


32| . and may happen because we begin with the system far from equilibrium. We 


integrated the classical equations as far as t = 60 and the solutions continued to be highly 
irregular. However, further investigation is beyond the scope of the present work. 

In Fig. in] we show the non-interacting ^ij values with an initial Fock state, averaged over 
2.88 X 10^ trajectories. We see that all values are positive semi-dehnite and irregular in time. 
The same correlations with x = 10“^ are presented in Fig. [TOl These were averaged over 
5.54 X 10® trajectories, and show clearly how interactions degrade the entanglement that can 
be detected by the Hillery-Zubairy measure. As expected, we see that this measure shows 
the onset of entanglement between the hrst well and the others at times which increase for 
the position along the chain. The entanglement between wells one and four with interactions 
is also seen to persist at later times than that between one and three, at least over this time 
interval. When we begin with an initial coherent state in the hrst well, only (^12 takes on a 
positive value over the evolution time, with ^12 ~ 1 around Jt = 6.35. 

There are six different bipartite EPR steering correlations between the hrst well and the 
other three. In the non-interacting case with an initial Fock state, they all take on positive 
values at some times, as shown in Fig. [11] and Fig. [12] As found with the two-and three-well 
systems, and Sjj were never found to be positive at the same time. This hints that 
asymmetric EPR steering may not be at all uncommon in nature, although it is not possible 
to make a dehnitive claim using only one type of measurement. 

The results for the EPR steering correlations for the interacting system with an initial 
Fock state are presented in Fig. [13] Out of the six possible bipartite correlations, only the 
three shown are ever signihcantly positive. S 41 momentarily attained a value of less than 
0.05 around Jt = 0.7, which would be unlikely to be of practical importance. With an initial 
coherent state, these EPR steering measures did not take on positive values over the time 
interval investigated. This demonstrates that the interactions mainly serve to degrade the 
correlations and that, if we wish to see entanglement and EPR steering in a Bose-Hubbard 
system, it is preferable to start with a non-classical initial state and minimise collisional 


interactions, possibly using Feshbach resonance techniques 


4l| . In these systems, a quantum 


initial state plays a greater role than in nonlinear optical systems, where, for example, a 
squeezed pump was predicted to improve quantum correlations in travelling-wave second 
harmonic generation The essential difference is that second harmonic generation is 


18 













FIG. 9: (Color online) The numerical non-interacting solutions for the Hillery-Zubairy criteria 
of the four-well model as a function of time, for J = 1, and A^i(O) = 200 in a Fock state, with 
-^ 2 ( 0 ) = A^3(0) = A^4(0) = 0. These solutions were averaged over 2.88 x 10^ stochastic realisations 
of the positive-P equations. 

perfectly capable of producing squeezed outputs even beginning from a completely classical 
pump. 

VII. CONCLUSIONS 

We have investigated the short time dynamics of the populations, entanglement and 
EPR steering in small Bose-Hubbard chains with a mesoscopic number of atoms initially 
inhabiting only the hrst well. The tunnelling is turned on at t = 0, which is equivalent to 
a sudden change of parameters. Using correlations that are suitable for number conserving 
processes, we hnd that bipartite entanglement and EPR steering spreads along the chain 
and gives a periodic signal in the two- and three-well cases without interatomic interactions, 
but only when the initial state in the hrst well is sub-Poissonian. In the cases where we 
could obtain analytical solutions, we saw that the existence of the entanglement and EPR 
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FIG. 10: (Color online) The numerical interacting solutions for the Hillery-Zubairy criteria of the 
four-well model as a function of time, for J = 1, y = 10“^ and A^i(O) = 200 in a Fock state, with 
-^ 2 ( 0 ) = ^^ 3 ( 0 ) = A^ 4 ( 0 ) = 0. These solutions were averaged over 5.54 x 10® stochastic realisations 
of the positive-P equations. 

steering signals depended completely on the number statistics in the first well being sub- 
Poissonian. The Fock states we used are maximally sub-Poissonian, so that the signals 
were optimised for these cases. When we began with an initial coherent state, we saw no 
evidence of either entanglement or EPR steering in these cases. This is not unexpected 
since the linear tunnelling is somewhat analogous to a time dependent beamsplitter and it 
is known that this cannot produce entangled outputs from classical inputs. In the four-well 
system, the solutions for the populations and the different correlations were irregular and 
non-periodic, even without interactions. In this case an initial Fock state resulted in strong 
signals for entanglement and EPR steering between the two end wells. 

The solutions for the well populations did not change signihcantly over the time scales 
investigated, being virtually indistinguishable in the two- and three-well systems. The four- 
well population solutions did begin to diverge from the non-interacting populations after 
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FIG. 11: (Color online) The numerical non-interacting solutions for the EPR steering criteria of 
the four-well model as a function of time, for J = 1 and Ni{0) = 200 in a Fock state, with 
-^ 2 ( 0 ) = A^3(0) = A^4(0) = 0. These solutions were averaged over 2.88 x 10® stochastic realisations 
of the positive-P equations. We again see that the system is asymmetric with regard to these 
measures of EPR steering. 

some time, but followed the same general form. The correlations for an initial Fock state 
however, were seriously affected by the interactions, only following the non-interacting so¬ 
lutions at short times. Some individual correlations which had signalled entanglement or 
EPR steering in the non-interacting case failed to do so once interactions were added. On 
the other hand, some correlations which failed to attain positive values for initial coherent 
states did become positive with interactions, although only over small time intervals and 
with small maximum values. In every case, the EPR steering correlations between any two 
wells were found to be different under a change of indices, with Sjj and Sj* never having 
positive values simultaneously. We cannot say if this would be the case for all possible EPR 
steering measurements, but any steering is asymmetric for the measure we used here. 

Overall, the initial quantum state in the hrst well was found to be far more important 
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FIG. 12: (Color online) The numerical non-interacting solutions for the EPR steering criteria S 14 
and S 41 , for J = 1, X = 10“^ and A^i(O) = 200 in a Fock state, with A^ 2 ( 0 ) = N^iO) = ^ 4 ( 0 ) = 0. 
These solutions were averaged over 2.88 x 10® stochastic realisations of the positive-P equations. 

in achieving the bipartite correlations than the interactions. This suggests that if Bose- 
Hubbard dynamics are to be used for quantum information purposes, it would be preferable 
to lower the collisional interactions as far as possible and concentrate on producing non- 
classical initial states. As the natural state of the atoms in an isolated well is a Fock state, 
this may provide signihcant advantages. 

Acknowledgments 

This research was supported by the Australian Research Council under the Future Fel¬ 
lowships Program (Grant ID: FT100100515). The author thanks Simon Haine for invaluable 


22 














FIG. 13: (Color online) The numerical interacting solutions for the EPR steering criteria of the 
four-well model as a function of time, for J = 1, x = 10“^ and A^i(O) = 200 in a Fock state, with 
-^ 2 ( 0 ) = ^^ 3 ( 0 ) = A^ 4 ( 0 ) = 0. These solutions were averaged over 5.54 x 10^ stochastic realisations 
of the positive-P equations. Apart from S 41 , which momentarily attained a value of less than 0.05 
around Jt = 0.7, those not shown remained negative. 

help with Mathematica. 
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